GenoTools: An Open-Source Python Package for Efficient Genotype Data Quality Control and Analysis.

GenoTools, a Python package, streamlines population genetics research by integrating ancestry estimation, quality control (QC), and genome-wide association studies (GWAS) capabilities into efficient pipelines. By tracking samples, variants, and quality-specific measures throughout fully customizable pipelines, users can easily manage genetics data for large and small studies. GenoTools' "Ancestry" module renders highly accurate predictions, allowing for high-quality ancestry-specific studies, and enables custom ancestry model training and serialization, specified to the user's genotyping or sequencing platform. As the genotype processing engine that powers several large initiatives including the NIH's Center for Alzheimer's and Related Dementias (CARD) and the Global Parkinson's Genetics Program (GP2). GenoTools was used to process and analyze the UK Biobank and major Alzheimer's Disease (AD) and Parkinson's Disease (PD) datasets with over 400,000 genotypes from arrays and 5000 sequences and has led to novel discoveries in diverse populations. It has provided replicable ancestry predictions, implemented rigorous QC, and conducted genetic ancestry-specific GWAS to identify systematic errors or biases through a single command. GenoTools is a customizable tool that enables users to efficiently analyze and scale genotype data with reproducible and scalable ancestry, QC, and GWAS pipelines.

Analysis of rare Parkinson's disease variants in millions of people.

Although many rare variants have been reportedly associated with Parkinson's disease (PD), many have not been replicated or have failed to replicate. Here, we conduct a large-scale replication of rare PD variants. We assessed a total of 27,590 PD cases, 6701 PD proxies, and 3,106,080 controls from three data sets: 23andMe, Inc., UK Biobank, and AMP-PD. Based on well-known PD genes, 834 variants of interest were selected from the ClinVar annotated 23andMe dataset. We performed a meta-analysis using summary statistics of all three studies. The meta-analysis resulted in five significant variants after Bonferroni correction, including variants in GBA1 and LRRK2. Another eight variants are strong candidate variants for their association with PD. Here, we provide the largest rare variant meta-analysis to date, providing information on confirmed and newly identified variants for their association with PD using several large databases. Additionally we also show the complexities of studying rare variants in large-scale cohorts.

Proxy-analysis of the genetics of cognitive decline in Parkinson's disease through polygenic scores.

Cognitive decline is common in Parkinson's disease (PD) and its genetic risk factors are not well known to date, besides variants in the GBA and APOE genes. However, variation in complex traits is caused by numerous variants and is usually studied with genome-wide association studies (GWAS), requiring a large sample size, which is difficult to achieve for outcome measures in PD. Taking an alternative approach, we computed 100 polygenic scores (PGS) related to cognitive, dementia, stroke, and brain anatomical phenotypes and investigated their association with cognitive decline in six longitudinal cohorts. The analysis was adjusted for age, sex, genetic ancestry, follow-up duration, GBA and APOE status. Then, we meta-analyzed five of these cohorts, comprising a total of 1702 PD participants with 6156 visits, using the Montreal Cognitive Assessment as a cognitive outcome measure. After correction for multiple comparisons, we found four PGS significantly associated with cognitive decline: intelligence (p = 5.26e-13), cognitive performance (p = 1.46e-12), educational attainment (p = 8.52e-10), and reasoning (p = 3.58e-5). Survival analyses highlighted an offset of several years between the first and last quartiles of PGS, with significant differences for the PGS of cognitive performance (5 years) and educational attainment (7 years). In conclusion, we found four PGS associated with cognitive decline in PD, all associated with general cognitive phenotypes. This study highlights the common genetic factors between cognitive decline in PD and the general population, and the importance of the participant's cognitive reserve for cognitive outcome in PD.

omicSynth: An open multi-omic community resource for identifying druggable targets across neurodegenerative diseases.

Treatments for neurodegenerative disorders remain rare, but recent FDA approvals, such as lecanemab and aducanumab for Alzheimer disease (MIM: 607822), highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use summary-data-based Mendelian randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer disease, 3 amyotrophic lateral sclerosis (MIM: 105400), 5 Lewy body dementia (MIM: 127750), 46 Parkinson disease (MIM: 605909), and 9 progressive supranuclear palsy (MIM: 601104) target genes passing multiple test corrections (pSMR_multi < 2.95 × 10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics, classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these, 69.8% are expressed in the disease-relevant cell type from single-nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development, and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as riluzole in Alzheimer disease. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community.

Multi-ancestry genome-wide association meta-analysis of Parkinson's disease.

Although over 90 independent risk variants have been identified for Parkinson's disease using genome-wide association studies, most studies have been performed in just one population at a time. Here we performed a large-scale multi-ancestry meta-analysis of Parkinson's disease with 49,049 cases, 18,785 proxy cases and 2,458,063 controls including individuals of European, East Asian, Latin American and African ancestry. In a meta-analysis, we identified 78 independent genome-wide significant loci, including 12 potentially novel loci (MTF2, PIK3CA, ADD1, SYBU, IRS2, USP8, PIGL, FASN, MYLK2, USP25, EP300 and PPP6R2) and fine-mapped 6 putative causal variants at 6 known PD loci. By combining our results with publicly available eQTL data, we identified 25 putative risk genes in these novel loci whose expression is associated with PD risk. This work lays the groundwork for future efforts aimed at identifying PD loci in non-European populations.

Sleep disturbances as risk factors for neurodegeneration later in life.

The relationship between sleep disorders and neurodegeneration is complex and multi-faceted. Using over one million electronic health records (EHRs) from Wales, UK, and Finland, we mined biobank data to identify the relationships between sleep disorders and the subsequent manifestation of neurodegenerative diseases (NDDs) later in life. We then examined how these sleep disorders' severity impacts neurodegeneration risk. Additionally, we investigated how sleep attributed risk may compensate for the lack of genetic risk factors (i.e. a lower polygenic risk score) in NDD manifestation. We found that sleep disorders such as sleep apnea were associated with the risk of Alzheimer's disease (AD), amyotrophic lateral sclerosis, dementia, Parkinson's disease (PD), and vascular dementia in three national scale biobanks, with hazard ratios (HRs) ranging from 1.31 for PD to 5.11 for dementia. These sleep disorders imparted significant risk up to 15 years before the onset of an NDD. Cumulative number of sleep disorders in the EHRs were associated with a higher risk of neurodegeneration for dementia and vascular dementia. Sleep related risk factors were independent of genetic risk for Alzheimer's and Parkinson's, potentially compensating for low genetic risk in overall disease etiology. There is a significant multiplicative interaction regarding the combined risk of sleep disorders and Parkinson's disease. Poor sleep hygiene and sleep apnea are relatively modifiable risk factors with several treatment options, including CPAP and surgery, that could potentially reduce the risk of neurodegeneration. This is particularly interesting in how sleep related risk factors are significantly and independently enriched in manifesting NDD patients with low levels of genetic risk factors for these diseases.

NeuroBooster Array: A Genome-Wide Genotyping Platform to Study Neurological Disorders Across Diverse Populations.

Genome-wide genotyping platforms have the capacity to capture genetic variation across different populations, but there have been disparities in the representation of population-dependent genetic diversity. The motivation for pursuing this endeavor was to create a comprehensive genome-wide array capable of encompassing a wide range of neuro-specific content for the Global Parkinson's Genetics Program (GP2) and the Center for Alzheimer's and Related Dementias (CARD). CARD aims to increase diversity in genetic studies, using this array as a tool to foster inclusivity. GP2 is the first supported resource project of the Aligning Science Across Parkinson's (ASAP) initiative that aims to support a collaborative global effort aimed at significantly accelerating the discovery of genetic factors contributing to Parkinson's disease and atypical parkinsonism by generating genome-wide data for over 200,000 individuals in a multi-ancestry context. Here, we present the Illumina NeuroBooster array (NBA), a novel, high-throughput and cost-effective custom-designed content platform to screen for genetic variation in neurological disorders across diverse populations. The NBA contains a backbone of 1,914,934 variants (Infinium Global Diversity Array) complemented with custom content of 95,273 variants implicated in over 70 neurological conditions or traits with potential neurological complications. Furthermore, the platform includes over 10,000 tagging variants to facilitate imputation and analyses of neurodegenerative disease-related GWAS loci across diverse populations. The NBA can identify low frequency variants and accurately impute over 15 million common variants from the latest release of the TOPMed Imputation Server as of August 2023 (reference of over 300 million variants and 90,000 participants). We envisage this valuable tool will standardize genetic studies in neurological disorders across different ancestral groups, allowing researchers to perform genetic research inclusively and at a global scale.

Long-Read Sequencing Resolves a Complex Structural Variant in PRKN Parkinson's Disease.

BACKGROUND: Parkin RBR E3 ubiquitin-protein ligase (PRKN) mutations are the most common cause of young onset and autosomal recessive Parkinson's disease (PD). PRKN is located in FRA6E, which is one of the common fragile sites in the human genome, making this region prone to structural variants. However, complex structural variants such as inversions of PRKN are seldom reported, suggesting that there are potentially unrevealed complex pathogenic PRKN structural variants. OBJECTIVES: To identify complex structural variants in PRKN using long-read sequencing. METHODS: We investigated the genetic cause of monozygotic twins presenting with a young onset dystonia-parkinsonism using targeted sequencing, whole exome sequencing, multiple ligation probe amplification, and long-read sequencing. We assessed the presence and frequency of complex inversions overlapping PRKN using whole-genome sequencing data of Accelerating Medicines Partnership Parkinson's disease (AMP-PD) and United Kingdom (UK)-Biobank datasets. RESULTS: Multiple ligation probe amplification identified a heterozygous exon three deletion in PRKN and long-read sequencing identified a large novel inversion spanning over 7 Mb, including a large part of the coding DNA sequence of PRKN. We could diagnose the affected subjects as compound heterozygous carriers of PRKN. We analyzed whole genome sequencing data of 43,538 participants of the UK-Biobank and 4941 participants of the AMP-PD datasets. Nine inversions in the UK-Biobank and two in AMP PD were identified and were considered potentially damaging and likely to affect PRKN expression. CONCLUSIONS: This is the first report describing a large 7 Mb inversion involving breakpoints outside of PRKN. This study highlights the importance of using long-read sequencing for structural variant analysis in unresolved young-onset PD cases. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Long-read sequencing resolves a complex structural variant in PRKN Parkinson's disease.

Background: PRKN mutations are the most common cause of young onset and autosomal recessive Parkinson's disease (PD). PRKN is located in FRA6E which is one of the common fragile sites in the human genome, making this region prone to structural variants. However, complex structural variants such as inversions of PRKN are seldom reported, suggesting that there are potentially unrevealed complex pathogenic PRKN structural variants. Objectives: To identify complex structural variants in PRKN using long-read sequencing. Methods: We investigated the genetic cause of monozygotic twins presenting with a young onset dystonia-parkinsonism using targeted sequencing, whole exome sequencing, multiple ligation probe amplification, and long-read. We assessed the presence and frequency of complex inversions overlapping PRKN using whole-genome sequencing data of AMP-PD and UK-Biobank datasets. Results: Multiple ligation probe amplification identified a heterozygous exon 3 deletion in PRKN and long-read sequencing identified a large novel inversion spanning over 7Mb, including a large part of the coding DNA sequence of PRKN. We could diagnose the affected subjects as compound heterozygous carriers of PRKN. We analyzed whole genome sequencing data of 43,538 participants of the UK-Biobank and 4,941 participants of the AMP-PD datasets. Nine inversions in the UK-Biobank and two in AMP PD were identified and were considered potentially damaging and likely to affect PRKN isoforms. Conclusions: This is the first report describing a large 7Mb inversion involving breakpoints outside of PRKN. This study highlights the importance of using long-read whole genome sequencing for structural variant analysis in unresolved young-onset PD cases.

Dopamine pathway and Parkinson's risk variants are associated with levodopa-induced dyskinesia.

Background: Levodopa-induced dyskinesia (LID) is a common adverse effect of levodopa, one of the main therapeutics used to treat the motor symptoms of Parkinson's disease (PD). Previous evidence suggests a connection between LID and a disruption of the dopaminergic system as well as genes implicated in PD, including GBA1 and LRRK2. Objectives: To investigate the effects of genetic variants on risk and time to LID. Methods: We performed a genome-wide association study (GWAS) and analyses focused on GBA1 and LRRK2 variants. We also calculated polygenic risk scores including risk variants for PD and variants in genes involved in the dopaminergic transmission pathway. To test the influence of genetics on LID risk we used logistic regression, and to examine its impact on time to LID we performed Cox regression including 1,612 PD patients with and 3,175 without LID. Results: We found that GBA1 variants were associated with LID risk (OR=1.65, 95% CI=1.21-2.26, p=0.0017) and LRRK2 variants with reduced time to LID onset (HR=1.42, 95% CI=1.09-1.84, p=0.0098). The fourth quartile of the PD PRS was associated with increased LID risk (ORfourth_quartile=1.27, 95% CI=1.03-1.56, p=0.0210). The third and fourth dopamine pathway PRS quartiles were associated with a reduced time to development of LID (HRthird_quartile=1.38, 95% CI=1.07-1.79, p=0.0128; HRfourth_quartile=1.38, 95% CI=1.06-1.78, p=0.0147). Conclusions: This study suggests that variants implicated in PD and in the dopaminergic transmission pathway play a role in the risk/time to develop LID. Further studies will be necessary to examine how these findings can inform clinical care.

Defining the causes of sporadic Parkinson's disease in the global Parkinson's genetics program (GP2).

The Global Parkinson's Genetics Program (GP2) will genotype over 150,000 participants from around the world, and integrate genetic and clinical data for use in large-scale analyses to dramatically expand our understanding of the genetic architecture of PD. This report details the workflow for cohort integration into the complex arm of GP2, and together with our outline of the monogenic hub in a companion paper, provides a generalizable blueprint for establishing large scale collaborative research consortia.

Polygenic risk score for Parkinson's disease and olfaction among middle-aged to older women.

INTRODUCTION: Olfactory impairment and Parkinson's disease (PD) may share common genetic and environmental risk factors. This study investigates the association of a PD polygenic risk score (PRS) with olfaction, and whether the associations are modified by environmental exposures of PM2.5, NO2, or smoking. METHODS: This analysis included 3358 women (aged 50-80) from the Sister Study with genetic data and results from the Brief Smell Identification Test (B-SIT) administered in 2018-2019. PD PRS was calculated using 90 single nucleotide polymorphisms. Olfactory impairment was defined with different B-SIT cutoffs, and PD diagnosis was adjudicated via expert review. We report odds ratios (ORs) and 95% confidence intervals (CIs) from multivariable logistic regression. RESULTS: As expected, PD PRS was strongly associated with the odds of having PD (OR highest vs. lowest quartile = 3.79 (1.64, 8.73)). The highest PRS quartile was also associated with olfactory impairment, with OR ranging from 1.24 (0.98, 1.56) for a B-SIT cutoff of 9 to 1.42 (1.04, 1.92) for a cutoff of 6. For individual B-SIT items, the highest PRS quartile was generally associated with lower odds of correctly identifying the odorant, albeit only statistically significant for pineapple (0.72 (0.56, 0.94), soap (0.76 (0.58, 0.99)) and rose (0.70 (0.54, 0.92)). The association of PD PRS with olfactory impairment was not modified by airborne environmental exposures or smoking. CONCLUSION: These preliminary data suggest that high PD genetic susceptibility is associated with olfactory impairment in middle-aged and older women.

Genome-Wide Meta-Analysis of Cerebrospinal Fluid Biomarkers in Alzheimer's Disease and Parkinson's Disease Cohorts.

BACKGROUND: Amyloid-ß, phosphorylated tau (p-tau), and total tau (t-tau) in cerebrospinal fluid are established biomarkers for Alzheimer's disease (AD). In other neurodegenerative diseases, such as Parkinson's disease (PD), these biomarkers have also been found to be altered, and the molecular mechanisms responsible for these alterations are still under investigation. Moreover, the interplay between these mechanisms and the diverse underlying disease states remains to be elucidated. OBJECTIVE: To investigate genetic contributions to the AD biomarkers and assess the commonality and heterogeneity of the associations per underlying disease status. METHODS: We conducted genome-wide association studies (GWASs) for the AD biomarkers on subjects from the Parkinson's Progression Markers Initiative, the Fox Investigation for New Discovery of Biomarkers, and the Alzheimer's Disease Neuroimaging Initiative, and meta-analyzed with the largest AD GWAS. We tested heterogeneity of associations of interest between different disease statuses (AD, PD, and control). RESULTS: We observed three GWAS signals: the APOE locus for amyloid-ß, the 3q28 locus between GEMC1 and OSTN for p-tau and t-tau, and the 7p22 locus (top hit: rs60871478, an intronic variant for DNAAF5, also known as HEATR2) for p-tau. The 7p22 locus is novel and colocalized with the brain DNAAF5 expression. Although no heterogeneity from underlying disease status was observed for the earlier GWAS signals, some disease risk loci suggested disease-specific associations with these biomarkers. CONCLUSIONS: Our study identified a novel association at the intronic region of DNAAF5 associated with increased levels of p-tau across all diseases. We also observed some disease-specific genetic associations with these biomarkers. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Genetic meta-analysis of levodopa induced dyskinesia in Parkinson's disease.

The genetic basis of levodopa-induced-dyskinesia (LiD) is poorly understood, and there have been few well-powered genome-wide studies. We performed a genome-wide survival meta-analyses to study the effect of genetic variation on the development of LiD in five separate longitudinal cohorts, and meta-analysed the results. We included 2784 PD patients, of whom 14.6% developed LiD. We found female sex (HR = 1.35, SE = 0.11, P = 0.007) and younger age at onset (HR = 1.8, SE = 0.14, P = 2 × 10-5) increased the probability of developing LiD. We identified three genetic loci significantly associated with time-to-LiD onset. rs72673189 on chromosome 1 (HR = 2.77, SE = 0.18, P = 1.53 × 10-8) located at the LRP8 locus, rs189093213 on chromosome 4 (HR = 3.06, SE = 0.19, P = 2.81 × 10-9) in the non-coding RNA LINC02353 locus, and rs180924818 on chromosome 16 (HR = 3.13, SE = 0.20, P = 6.27 × 10-9) in the XYLT1 locus. Based on a functional annotation analysis on chromosome 1, we determined that changes in DNAJB4 gene expression, close to LRP8, are an additional potential cause of increased susceptibility to LiD. Baseline anxiety status was significantly associated with LiD (OR = 1.14, SE = 0.03, P = 7.4 × 10-5). Finally, we performed a candidate variant analysis of previously reported loci, and found that genetic variability in ANKK1 (rs1800497, HR = 1.27, SE = 0.09, P = 8.89 × 10-3) and BDNF (rs6265, HR = 1.21, SE = 0.10, P = 4.95 × 10-2) loci were significantly associated with time to LiD in our large meta-analysis.

Genome-wide meta-analysis of CSF biomarkers in Alzheimer's disease and Parkinson's disease cohorts.

Background: Amyloid beta (Aß), phosphorylated tau (p-tau), and total tau (t-tau) in cerebrospinal fluid are established biomarkers for Alzheimer's disease (AD). In other neurodegenerative diseases, such as Parkinson's disease (PD), these biomarkers have also been found to be altered, and the molecular mechanisms responsible for these alterations are still under investigation. Moreover, the interplay between these mechanisms and the diverse underlying disease states remains to be elucidated. Objectives: To investigate genetic contributions to the AD biomarkers and assess the commonality and heterogeneity of the associations per underlying disease status. Methods: We conducted GWAS for the AD biomarkers on subjects from the Parkinson's Progression Markers Initiative (PPMI), the Fox Investigation for New Discovery of Biomarkers (BioFIND), and the Alzheimer's Disease Neuroimaging Initiative (ADNI) and meta-analyzed with the largest AD GWAS.[7] We tested heterogeneity of associations of interest between different disease statuses (AD, PD, and control). Results: We observed three GWAS signals: the APOE locus for Aß, the 3q28 locus between GEMC1 and OSTN for p-tau and t-tau, and the 7p22 locus (top hit: rs60871478, an intronic variant for DNAAF5 , also known as HEATR2 ) for p-tau. The 7p22 locus is novel and co-localized with the brain DNAAF5 expression. While no heterogeneity from underlying disease status was observed for the above GWAS signals, some disease risk loci suggested disease specific associations with these biomarkers. Conclusions: Our study identified a novel association at the intronic region of DNAAF5 associated with increased levels of p-tau across all diseases. We also observed some disease specific genetic associations with these biomarkers.

Application of Aligned-UMAP to longitudinal biomedical studies.

High-dimensional data analysis starts with projecting the data to low dimensions to visualize and understand the underlying data structure. Several methods have been developed for dimensionality reduction, but they are limited to cross-sectional datasets. The recently proposed Aligned-UMAP, an extension of the uniform manifold approximation and projection (UMAP) algorithm, can visualize high-dimensional longitudinal datasets. We demonstrated its utility for researchers to identify exciting patterns and trajectories within enormous datasets in biological sciences. We found that the algorithm parameters also play a crucial role and must be tuned carefully to utilize the algorithm's potential fully. We also discussed key points to remember and directions for future extensions of Aligned-UMAP. Further, we made our code open source to enhance the reproducibility and applicability of our work. We believe our benchmarking study becomes more important as more and more high-dimensional longitudinal data in biomedical research become available.

Genetic meta-analysis of levodopa induced dyskinesia in Parkinson's disease.

Importance: Forty percent of Parkinson's disease patients develop levodopa-induced-dyskinesia (LiD) within 4 years of starting levodopa. The genetic basis of LiD remains poorly understood, and there have been few well powered studies. Objective: To discover common genetic variants in the PD population that increase the probability of developing LiD. Design setting and Participants: We performed survival analyses to study the development of LiD in 5 separate longitudinal cohorts. We performed a meta-analysis to combine the results of genetic association from each study based on a fixed effects model weighting the effect sizes by the inverse of their standard error. The selection criteria was specific to each cohort. We studied individuals that were genotyped from each cohort and that passed our analysis specific inclusion criteria. Main Outcomes and Measures: We measured the time for PD patients on levodopa treatment to develop LiD as defined by reaching a score higher or equal than 2 from the MDS-UPDRS part IV, item 1, which is equivalent to a range of 26%-50% of the waking time with dyskinesia. We carried out a genome-wide analysis of the hazard ratio and the association of genome-wide SNPs with the probability of developing LiD using cox proportional hazard models (CPH). Results: This study included 2,784 PD patients of European ancestry, of whom 14.6% developed LiD. Consistent with previous studies, we found female gender (HR = 1.35, SE = 0.11, P = 0.007) and younger age at onset (HR = 1.8, SE = 0.14, P = 2 × 10 -5 ) to increase the probability of developing LiD. We identified three loci significantly associated with time-to-LiD onset. rs72673189 on chromosome 1 (HR = 2.77, SE = 0.18, P = 1.53 × 10 -8 ) located in the LRP8 locus, rs189093213 on chromosome 4 (HR = 3.06,, SE = 0.19, P = 2.81 × 10 -9 ) in the non-coding RNA LINC02353 locus, and rs180924818 on chromosome 16 (HR = 3.13, SE = 0.20, P = 6.27 × 10 -9 ) in the XYLT1 locus. Subsequent colocalization analyses on chromosome 1 identified DNAJB4 as a candidate gene associated with LiD through a change in gene expression. We computed a PRS based on our GWAS meta-analysis and found high accuracy to stratify between PD-LID and PD (AUC 83.9). We also performed a stepwise regression analysis for baseline features selection associated with LiD status. We found baseline anxiety status to be significantly associated with LiD (OR = 1.14, SE = 0.03, P = 7.4 × 10 -5 ). Finally, we performed a candidate variant analysis and found that genetic variability in ANKK1 ( rs1800497 , Beta = 0.24, SE = 0.09, P = 8.89 × 10 -3 ) and BDNF ( rs6265 , Beta = 0.19, SE = 0.10, P = 4.95 × 10 -2 ) loci were significantly associated with time to LiD in our large meta-analysis. Conclusion: In this association study, we have found three novel genetic variants associated with LiD, as well as confirming reports that variability in ANKK1 and BDNF loci were significantly associated with LiD probability. A PRS nominated from our time-to-LiD meta-analysis significantly differentiated between PD-LiD and PD. In addition, we have found female gender, young PD onset and anxiety to be significantly associated with LiD.

Large-scale rare variant burden testing in Parkinson's disease.

Parkinson's disease has a large heritable component and genome-wide association studies have identified over 90 variants with disease-associated common variants, providing deeper insights into the disease biology. However, there have not been large-scale rare variant analyses for Parkinson's disease. To address this gap, we investigated the rare genetic component of Parkinson's disease at minor allele frequencies <1%, using whole genome and whole exome sequencing data from 7184 Parkinson's disease cases, 6701 proxy cases and 51 650 healthy controls from the Accelerating Medicines Partnership Parkinson's disease (AMP-PD) initiative, the National Institutes of Health, the UK Biobank and Genentech. We performed burden tests meta-analyses on small indels and single nucleotide protein-altering variants, prioritized based on their predicted functional impact. Our work identified several genes reaching exome-wide significance. Two of these genes, GBA1 and LRRK2, have variants that have been previously implicated as risk factors for Parkinson's disease, with some variants in LRRK2 resulting in monogenic forms of the disease. We identify potential novel risk associations for variants in B3GNT3, AUNIP, ADH5, TUBA1B, OR1G1, CAPN10 and TREML1 but were unable to replicate the observed associations across independent datasets. Of these, B3GNT3 and TREML1 could provide new evidence for the role of neuroinflammation in Parkinson's disease. To date, this is the largest analysis of rare genetic variants in Parkinson's disease.

Genetic risk factor clustering within and across neurodegenerative diseases.

Overlapping symptoms and co-pathologies are common in closely related neurodegenerative diseases (NDDs). Investigating genetic risk variants across these NDDs can give further insight into disease manifestations. In this study we have leveraged genome-wide single nucleotide polymorphisms and genome-wide association study summary statistics to cluster patients based on their genetic status across identified risk variants for five NDDs (Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, Lewy body dementia and frontotemporal dementia). The multi-disease and disease-specific clustering results presented here provide evidence that NDDs have more overlapping genetic aetiology than previously expected and how neurodegeneration should be viewed as a spectrum of symptomology. These clustering analyses also show potential subsets of patients with these diseases that are significantly depleted for any known common genetic risk factors suggesting environmental or other factors at work. Establishing that NDDs with overlapping pathologies share genetic risk loci, future research into how these variants might have different effects on downstream protein expression, pathology and NDD manifestation in general is important for refining and treating NDDs.